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ABSTRACT 

We propose a new method to locate the 7-ray-emitting positions i? 7 from the measured 
time lags r Q b of 7-ray emission relative to broad emission lines. The method is also 
applicable to lower frequencies. i? 7 depends on parameters r b, -Rblr, and 9, where 
-Rblr is the size of broad-line region, Vd is the travelling speed of disturbances down 
the jet and 9 is the viewing angle of the jet axis to the line of sight. As r Q b = 0, 
T b < or r b > 0, the broad lines zero-lag, lag or lead the 7-rays, respectively. It 
is applied to 3C 273, in which the lines and the radio emission have enough data, 
but the 7-rays have not. We find r Q b < and r Q b > for the 5, 8, 15, 22 and 37 
GHz emission relative to the broad lines Ha, H/3 and H7. The lag may be positive or 
negative, however current data do not allow to discriminate between the two cases. 
The measured lags are on the order of years. For a given line, Tbb generally decreases 
as radio frequency increases. This trend most likely results from the radiative cooling 
of relativistic electrons. The negative lags have an average of T b = —2.86 years for the 
37 GHz emission, which represents that the lines lag the radio emission. The positive 
lags have r b = 3.20 years, which represents that the lines lead the radio emission. 
We obtain the radio emitting positions i? ra dio = 0.40-2.62 pc and i? ra dio = 9.43-62.31 
pc fo r the negative and posi tive la gs, respectively . From the constraint of i? 7 < i? ra dio 
(e.g. iDermer fc Schlickeiserl Il994 Uorstad et al.l l200lD . we have i? 7 < 0.40-2.62 pc 
for the negative lags. For the positive lags, 4.67-30.81 < i? 7 < 9.43-62.31 pc. These 
estimated i? 7 are consistent with those of other researches. These agreements confirm 
the reliability of the method and assumptions. The method may be also applicable to 
BL Lacertae objects, in which broad lines were detected. 

Key words: 7-rays: theory - galaxies: active - galaxies: jets - quasars: emission lines 
- quasars: individual: 3C 273. 



1 INTRODUCTION 

The 7-rays of blazars are generally believed to be generated 
by inverse Compton (IC) emission from a relativistic jet ori- 
ented at a small angle to the line of sight (|Blandford fc Reea 
1 19781 ). Two of the most accepted scenarios for broad-band 
emission from radio to 7-rays are the synchrotron self- 
Co mpton (SSC) and exte rnal Compton (EC) models (see 
e.g. IChisellini et al.l 1 19981 ). The broad-band spectral en- 
ergy distributio ns (SEDs) of blaza r s consist of two broa d 
bumps (see e.g. iFossati et a"Hll998l ; iGmsellini et aflll998ft . 
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The first component is from the synchrotron process, and 
the second one, generally peaking at the 7-ray regime, is 
generated by the IC emission of the same electron pop- 
ulation responsible for the synchrotr on emission (see e.g. 
iGhisellini etafl ll99Sl ; iBottcherl 1 19991 ). However, the posi- 
tions of 7-ray-emitting regions are still an open and contro- 
versial issue in the researches on blazars. It was suggested 
that 7-rays are produced within broad-line region (BLR) 
and that the 7-ray-emitting pos itions R y range roughly be - 
tween 0.03 and 0.3 pa r sec jpc) (|Ghisellini fc Madau|[l99l ). 
iBlandford fc Levinsonl l|l995l ) also suggested a sub-pc 7-ray- 
emitting region. It was argued that the radiative plasma 
in relativistic jets of powerful blazars are i nside the BLR 
(|Georganopoulos. Kirk fc Mastichiadisl l200ll ) . On the con- 
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trary, it was also argue d that the 7-ray-emitting regions 
are outside the BLR (|Lindfors. Valtaoia fe Tiirlerl 120051 ; 
ISokolov fc Marscher|[2005n ~ Internal absorp tion for 10 GeV- 
1 TeV 7-rays were used to constrain Rj |Liu fc Bail 120061 : 
iLiu. Bai fc Mal2008l ; [Bai. Liu fc Mall2009l ). Variability of the 
high energy flux indicates that the 7-ray-emitting positions 
cannot be too distant from the ce nt ral supermassive black 
hole (jGhisellini fc Madaul Il996l ; IGhisellini fc Tavecchiol 
2009), while the photon-photon absorption implies that the 
emitting positions cannot be too close to the black hole and 
its accretion di s c dGhisellini fc Madau 



2008 



1996; IWa 



Liu fc Baill2006l: ILiu .Bai fc Mall2008l: Isitarek fc Bednare 



Bai, Liu fc Mai 120091; IGhisellini fc Tavecchiol I200S 



Tavecchio fc Mazinll2009i ). Bracketed by the two limits, one 



obtains a few hundreds of Schwarzschild radii as the pre- 
ferred jet location where m os t of the dissipation occurs 
(Ghisellini & Tavecchiol 120091 ). IJorstad et all |200lD estab- 
lished a connection between ejections of superluminal ra- 
dio knots and 7-ray outbursts observed by EGRET. They 
concluded that the radio and 7-ray events are orig i natin 
from the same region of a relativistic jet. [Wagnerl (|200 
found that the blazar 7-ray emission might depend on the 
mass of the central black hole and that very high energy 
(VHE) 7-ray-emitting active galactic nuclei (AGN s) have 
the bl ack hole masses larger than 10 8 solar masses. iBlooml 

(2008) c onfirmed the radi o and 7-ray correlation of EGRET 
blazars. iThompsonl (2008) reported that a central feature of 
the EGRET results is the high degree of variability seen in 
many 7-ray sources, indicative of the powerful central en- 
gines at work in objects visible to 7-ray telescopes. 

The operation of Fermi Gamma Ray Spa ce Telescope 
(_Ferrm) -La rge Area Telescope (LAT) (see e.g. lAbdo et al.l 
I2009l . l2010al l presents an exceptional opportunity for under- 
standing the central engines operating in blazars. The First 
LAT AGN Catalog includes 709 AGNs, comprising 300 BL 
Lacertae objects (BL Lacs), 296 flat spectrum radio quasars 
(FSR Qs), 41 AGNs of oth er types and 72 AGNs of unknown 
types (Abdo et al.l l2010al ). The distribution of 7-ray pho- 
ton spect ral index is found to co rrelate strongly with blazar 
subclass (|Abdo et alj l2010al lc]). lAbdo et"aH (|2010bl ) found 
for Fermi blazars that variation amplitudes are larger for 
FSRQs a nd low/intermediate synchrotron frequen cy peaked 
BL Lacs. IGhisellini. Tavecchio fc Ghirlandal (|20Qgh obtained 
a few hundreds of Schwarzschild r a dii as the preferred jet 
dissipation region. IGhisellini et~ai1 (|2010l ) found for bright 
Fermi blazars t hat the jet dissipatio n region is within the 
BLR for FSRQs. IKovalev et alj (|2009l ) identified the pc-scale 
radio core as a like l y loca tion for both the 7-ray and radio 
flares. ISikora et aD £2009) suggested that the blazar emis- 
sion zone is located at pc-scale distances from the nucleus, 
and they also proposed that the pc-scale blazar activity can 
be occasionally accompa nied by dissipative events taking 
place at sub-pc distances. IGhisellini. Maraschi fc Tavecchiol 

(2009) argued that the bulk of the most luminous blazars 
already detected by Fermi should be charact erized by large 
black hole masses, around 10 9 solar masses. IKovalev et al.l 
(2009) found for Fermi blazars that the 7-ray photon flux 
correlates with the compact radio density flux. 

Disturbances in the central engine are likely trans- 
ported down the relativistic jets. This was supported by 
observations that dips in the X-ray emission are followed 
by ejections of bright superluminal knots in the radio 



jets of AGNs and mi c roquasars (e.g. Marschcr et al.l |2002| ; 
IChateriee et all 120091 ; lArshakian et all l2010l )~ The events 
in the central engine, where the X-rays are produced, will 
have a direct effect o n the events in the r adio jets (e.g. 
iMarscher et~all |2002| ; IChateriee et all 120091 ). The disc-jet 
connection was suggested by correlations of emission line lu- 
minosity and radio power of je t s for various samples of AGNs 
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1995; Cclotti, Padovani & Ghisellini 1997; Serjeant et al 



19981; ICao fc Jiand . 
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1999. 



20011 ; IWang fc Hoi 120031 ; 

The relativistic jets can be 
ejected from inner accre tion disc in the vicinity of the cen- 
tral b l ack hole (see e.g. Penrose! 19691: Blandford fc Znaiek 



1977i IBlandford fc Payne! 1 19821 ; IMeier. Koide fc Uchida 
20011 ). The close connection between accretion disc and 



jets indicates that the central disturbances are likely trans- 
ported down the jets. For example, the central disturbances 
that drive variations of the central ionizing continuum 
may be transported with the ou t ward Alfven waves (see 
e.g. IMeier. Koide fc Uchidal l200ll ; iKoide et all |2002|). By 
magn etohydrodynamic mechanisms (|Kudoh fc Shibatal 
1 1999f ) . the transported disturbances may be converted into 
the local energies of plasma far from the central black hole. 
Thus the disturbances in the central engine would influence 
the 7-rays emitted by the relativistic jet aligned with the 
line of sight. At the same time, these disturbances can lead 
to the variations of the central ionizing continuum that 
drives broad emission lines from BLR. Hence, the broad 
lines from the BLR and the 7-rays from the relativistic jet 
may be both coupled to the disturbances in the central 
engine. 

Based on the photoionization assumption and the 
time lags between broad lines and continuum, rever- 
beration mapping observations are able to determine 
the sizes of the BL Rs for type 1 AGNs (see e.g . 



Kaspi fc Netzen 1999 ; Wandel, Peterson fc MalkarJ 1999; 
Kaspi et all |2000|. 120051. 120071 : iPeterson et all l 200oi . \2QQ4 . 



20051 ; IVestergaard fc Peterson! I2006IT According to the re - 
verberation mapping model (e.g. IBlandford fc McKeelll982l ). 
the variations of broad lines can reflect the disturbances in 
the central engines of blazars, even though the beamed emis- 
sion from the relativistic jet strongly affects the real ionizing 
continuum from accretion disc so that no appropriate con- 
tinuum could be used as the reference to estimate the time 
lags relative to the broad lines. If the broad lines from the 
BLR and the 7-rays from the relativistic jet are both coupled 
to the disturbances in the central engine, the disturbances 
should similarly influence both variations of the 7-rays and 
the broad lines. Thus there should be correlations and time 
lags between the broad lines and the 7-rays. The time lags 
should be related to R 1 . In the paper, we attempt to locate 
R 1 from the time lags between both variations of the broad 
lines from the BLR and the 7-rays from the relativistic jet. 



2 METHOD 

The broad lines from the BLR and the 7-rays from the 
relativistic jet may be both coupled to the disturbances 
in the central engine. Thus the disturbances could sim- 
ilarly influence both variations of the broad lines and 
the 7-rays emitted by the jet aligned with the line of 
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sight. The outbursts seen in light curves are physically 
linked to the ejections of superlu minal radio knots (e.g. 
ITiirler. Courvoisier fc Paltani 2000). The events in the cen- 
tral engine wi ll have a direct effect on the events in th e 
radio jets (e.g. iMarscher et alJliool IChateriee et al.ll2009h . 
Thus it is likely that the 7-ray outbursts are caused by the 
disturbances from the central engine, but not the local dis- 
turbances produced in the jets. Hence, there should be cor- 
relations between the 7-ray outbursts and the variations of 
broad lines. It is possible that there are time lags in the cor- 
relations and the time lags are related to the positions i? 7 . 
Thus R 7 could be located by the time lags between the 7-ray 
outbursts and the variations of broad lines. In the future, the 
quasi-simultaneous observations of 7-rays with Fermi/LAT 
and broad lines with optical telescopes on the order of years 
may be employed to test this expectation. The method of 
locating _R 7 is also applicable to infrared, optical and radio 
emission. 

First, it is assumed that the disturbances in the central 
engine are transported outward by some process, and the 
disturbances could similarly influence both variations of the 
broad lines from the BLR and the 7-rays from the relativis- 
tic jet aligned with the line of sight. Second, we assume a 
simple geometry (see Fig. la) that is similar to the schemes 
in the c lassical reverberation mapping of the broad lines 
(see e.g. Kaspi & Nctzcr 1999; Wandcl, Peterson & Malkan 
19991 ; iKaspi et all |2000| . |2005| . 120071 ; iPeterson et~ai1 |2000| . 



2004 I2005I ; IVestergaard fc Petersorfeod l. The differences 
between this method and the classical mapping are as fol- 
lows. In the classical mapping, a part of ionizing continuum 
drives the broad lines from the BLR, and another directly 
reaches observers. Thus the broad lines lag the continuum 
due to light travelling time effects, and the time lag corre- 
sponds to the size of the BLR. In the method proposed here, 
the ionizing photon signals directly detected by telescopes 
in the classical mapping are replaced with the transporting 
disturbances from the central engine down the jet and then 
with the jet emission signals at _R 7 . Thus the method would 
be valid even if the BLR has a complex configuration, such 
as a spherical shell. For simplicity for graphic illustration, 
we choose the BLR to be a ring (see Fig. la). The trans- 
porting speed of disturbances down the jet cannot exceed 
the speed of light. 

Because R 7 is unknown, the 7-rays may lag, lead or 
zero-lag the broad lines. First, we try to find the position 
where the 7-rays are produced and zero-lag the broad lines. 
As the disturbances reach point G (see Fig. la), where the 
7-rays are produced, i.e. _R 7 = AG, the ionizing continuum 
photons travel from point A to B and the line photons travel 
from point B to / in time interval R~,/vd- In the case of zero- 
lag (hereafter Case A), we have (Rblr + R 7 )/c = R 7 /vd, 
and then 



R~, 



Rblr ■ 



Vd 



(1) 



where Vd is the travelling speed of disturbances down the jet 
and c is the speed of light. For Case A, R 7 > Rblr. 

Within segment AG, the lines will lag the 7-rays. As the 
disturbances reach point F, where the 7-rays are produced, 
i.e. i? 7 = AF, the ionizing continuum photons reach point 
E. The light travelling time effects for the ionizing photons 
from point E to B and the line photons from B to K result in 



the time lag of the lines relative to the 7-rays (see Fig. la). In 
this case (hereafter Case B), we have R 7 + -Rblr — i? 7 c/ud = 
CT"ob/(l + z ), an d then 



-R 7 — 



RbLR. 



r ..h 
1+j 



(2) 



where z is the redshift of source, r b > and r b is the 
observed time lag of the broad lines relative to the 7-rays. 

Outside segment AG, the 7-rays will lag the lines. As 
the disturbances reach point H, where the 7-rays are pro- 
duced, i.e. R 7 = AH, the line photons reach point L. The 
light travelling time effects for the 7-ray photons from point 
J to L, i.e. from point H to M, result in the time lag of the 
7-rays relative to the lines (see Fig. la). In this case (here- 
after Case C), we have R 1 c/vd — R 7 + Rblr + cr b/(l + z), 
and then 



R-, 



R B LR + 



(3) 



where r b > and r b is the observed time lag of the 7-rays 
relative to the broad lines. Equations (1), (2) and (3) can 
be unified into 



R- 



Rblr + %k 
■f -1 : 



(4) 



where r D b is the observed time lags of the 7-rays relative to 
the broad lines and is zero, negative or positive. As r b = 
(Case A), equation (4) becomes equation (1). As r D b < 
(Case B), equation (4) becomes equation (2). As r D b > 
(Case C), equation (4) becomes equation (3). Once r b, 
-Rblr and Vd are known, R 7 can be obtained from equations 
(1)"(4). 

The calculations above are under the condition of 9 = 0, 
where 9 is the angle between the jet axis and the line of sight. 
In fact, the approaching relativistic jets of blazars are ori- 
ented at a small angle to the line of sig ht (|Blandford fc Reesi 
1 19T8T ) . Thus 9^0. Considering the actual inclination of the 
jet axis with respect to the line of sight (see Fig. lb), we re- 
deduce the expression of Case C. As the disturbances reach 
point H , where the 7-rays are produced, the line photons 
reach point L. The light travelling time effects for the 7-ray 
photons from point O to L, i.e. from point H to S, result in 
the time lag of the 7-rays relative to the lines (see Fig. 16). In 
this case, we have R^c/vd = -Rblr + -R 7 cos 9 + cr b/(l + 2), 
and then 



R~, 



Rblr + gg 
-^-cos(9 ' 



(•») 



It is obvious that equation (5) contains Cases A, B and C. As 
9 = 0, equation (5) becomes equation (4). The observed line 
photons are from the BLR of ring, and then the observed lag 
is an ensemble average over all points of the ring. For point 
T, we have Rjc/vd = -Rblr + TW + R 7 cos 8 + cr b/(l + z) 
and TW = -Rblr sin a sin 8 (see Fig. lb), and then 



R~, 



Rblr(1 + sin a sin 9) + 



- cos 9 



(6) 



where a is the angle between AB and AT and varies from 
to 2n. For a given source, i.e. given R 7 , -Rblr, Vd, 9 and 
z, r b varies with a, i.e. r b = r Q b(a). Calculating ensemble 
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Figure 1. Sketch of the geometry assumed, and it is similar to 
that used in the reverberation mapping method of broad emis- 
sion lines. -Rblr is the size of BLR. (a) the angle between the 
line of sight and the jet axis 8 = 0. (b) 8 ^ 0. The jet axis 
is perpendicular to the plane of BLR. (6) the planes ABU and 
HNOP are perpendicular to the line of sight, v TV _L AV, :. 
TV = AT sin a = .Rblr sin a. v TW \\ BL, :. TW _L WV and 
TW _L AV. TW _L AV and TV _L AV give WV _L AV. v 
TV _L AV and WV _L AV, :. ZTVW = 8. v ZTVW = 8 and 
ZTWV = w/2, :. TW = TV sin 8 = i? B LR.sinasm0. 

average over a in equation (6), we have 

-Rblr. + ^ 
= ^-co s0 > (7) 

where (r b) is the ensemble average of r b(a) and is the 
measured time lag between the 7-ray and line light curves 
(and (1 + sin a sin 8) — Jg^{l + sin a sin 6)da/2iv — 1). It is 
obvious that equation (7) contains Cases A, B and C. As 
= 0, equation (7) becomes equation (4). In the following 
sections, we will calculate R y based on equation (7). 



3 APPLICATION TO 3C 273 

3C 273 was first id entified as a quasar at redshift z = 0.158 
bv lSchmidtl d 19631). It is one of the best studied AGNs in all 
bands (see e . g. iLichti et all 1 19951; I von Montignv et aT1ll997l ; 
ITiirler et all Il999l ; ISoldi et al.ll2008r i. The database siteFI of 
3C 273 provides a series of about 70 light curves as well as 
some spectra. The jet of 3C 273 is one- sided, with no sign s 
of emission from the counterjet side (|Unwin et alj [l98a ). 

1 http://isdc.unige.ch/3c273/ 



The blue-bump of 3C 273 is thermal con tinuum emission 
from the inner accretion disc l|Shieldsll 19781 ). Fe Ka lines ob- 
served in 3C 273 were shown to be from an accretion disc 
around a supermassive black hole (lYaqoob fe Serlemitsosl 
2000; iTorres et alJ l~2003). If the assumptions in the method 
are correct, it is expected for 3C 273 that there should exist 
time lags between the broad lines from the BLR and the 
7-rays from the relativistic jet. 

3.1 Data of 3C 273 

This pape r makes use of the 3C 273 database hosted by 
the ISDC i|Tiirler et alj|l999h . This database is one of the 
most complete multi-wavelength databases currently avail- 
able for one AGN. For 3C 273, the 7-ray light curves are 
very sparsely sampled and/or the error bars of 7-ray fluxes 
are also very large. Of course with Fermi takin g data since 
June 2 008 this is no longer the case (see e.g. lAbdo et al.l 
l2010bl ). However, the sampling over the dates of interest, 
and, in general over the timescales of interest is still sparse 
(see the line light curves in the following paragraphs). Thus 
it should be unreliable to employ the 7-ray light curves 
to estimate the time lags. Hence for these blazars lacking 
adequate 7-ray light curves, the synchrotron emission, es- 
pecially the radio emission, could be used to derive the 
time lags relative to the broad lines. The synchrotron flares 
from the relativistic jet dominate energy output from ra- 
dio to millimeter and extend up to the infrared-optical 
regimes llRobson et al.l fl993l : ITiirler. Courvoisier fe Paltanil 
|2000| ; ISoldi et all 120081 ). Radio light curves are better than 
millimeter ones in the samplings and the features of syn- 
chrotron flares. Gamma-ray detections correspon d to rising 
radio fluxes (e.g. lUlrich. Maraschi fe Urrv||l"997l ). Thus the 
radio light curves are adopted to be analyzed. Though the 
synchrotron flares extend up to the infrared-optical regimes, 
the flares are sparse (denoted by red color in the 3C 273 
database) . In the intervals among these flares, the infrared- 
optical light curves may be contaminated by other emission 
components. Thus it is expected that there should be no sig- 
nificant features in the cross correlation functions between 
the broad lines and the light curves. Considering that the 
synchrotron emission peaks aroun d the infrared band for 3C 
273 (see e.g. iGhisellini et al.lll998T l , the infrared light curves 
are also adopted to be analyzed and to test the expecta- 
tion above. In the light curves adopted here, only good data 
(Flag>=0) are adopted. 

Light curves of 5, 8, 15, 22 and 37 GHz are taken from 
the 3C 273 databas e, for these light curves have enough data 
ijTiirler et al.|[l999l ). These radio light curves are presented 
in Fig. 2. There are four distinct outbursts in each radio 
light curve after ~ 1980 (see Fig. 2). The data of 22 and 
37 GHz are adopted from ~ 1980, for observations are very 
sharply sampled before ~ 1980 (see Fig. 2). From 1973 to 
1978, there is a gap of 5 years without observations of 5 
GHz. Considering the four distinct outbursts in each radio 
light curve, the data from ~ 1980 are adopted for the light 
curve of 5 GHz. Also, the data from ~ 1980 are adopted for 
the light curves of 8 and 15 GHz. The sampling rates of 5, 
8, 15, 22 and 37 GHz are 29, 40, 40, 44 and 46 times per 
year for the adopted data, respectively. 

Light curves of bro ad lines Ha, H/3 and H7 are taken 
from lKaspi et all l|2000l ). and all the data in the light curves 
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Figure 2. Light curves of the radio emission and the Balmer lines. The y-axis is in units of Jy. For lines, the y-axis is in units of 
10 -14 erg cm -2 s _1 . A, B, C and D denote the four distinct outbursts after ~ 1980. 



are adopted to calculate the time lags relative to the ra- 
dio emission. The sampling rates of the lines are around 5 
times per year. The line light curves are also presented in 
Fig. 2, wherein all the light curves show the same time in- 
terval. It can be seen in Fig. 2 that the line light curves are 
sharply sampled relative to the radio light curves adopted. 
Especially, there is only one data point in each of several 
valley-bottoms in the Ha and H/3 light curves, and the trend 
of variations can be changed if the points in these valley- 
bottoms are excluded from the Ha and H/3 light curves (see 
Fig. 2). The sharp samplings will affect the choice of analy- 
sis method for the cross correlation between the broad lines 
and the radio emission. The infrared light curves of J, H, 
K and L bands are also adopted from the 3C 273 database. 
The data before 1998 are adopted for the J, H and L light 
curves. The data before 2000 are adopted for the K light 
curve. 



3.2 Analysis of time lags 

Cross-correlation function (CCF) analysis is a standard 
technique in time series analysis for finding time lags be- 
tween light curves at different wavelengths, and the defi- 
nition of the CCF assumes that the light curves are uni- 
formly sampled. However, in most cases the sampling is not 
uniform. T he interpolated cros s corre lation function (ICCF) 
method of iGaskell fe Peterson! (|l987f ) uses a linear interpo- 
lation scheme to determine the missing data in the light 
curves. On the other hand, the discrete correlation function 
(DCF; lEdelson fe Krolik!ll988h can utilize a binning scheme 
to approximate the missing data. Apart from the ICCF and 



DCF, there is another method of estimating the CCF in the 
case of nonuniformly sampled light cu rves, the z-transf ormed 
discrete correlation function (ZDCF; lAlexanderH l997'). The 
ZDCF is a binning type of method as an improvement of 
the DCF technique, but it has a notable feature in that 
the data are binned by equal population rather than equal 
bin width as in the DCF. It has been shown in practice 
that the ZDCF is more robust than the ICCF and the DCF 
when ap plied to sparsely and unequally sampled l ight curves 

I 



Edelson et al 



y and unequally sampled light curves 
Il99d : iGiveon et all 1 19991 ; iRov et all 



I2000T ). iLiu et all IMPOST ) analyzed the ZDCFs between un- 
equally sampled light curves of AGNs, and they obtained 
inter-band time lags. In practice, the ZDCF is applicable 
and reliable to analyze the unequally sampled light curves. 
Thus the ZDCF will be calculated in the paper because of 
the sharp samplings of the Balmer lines. 

In general, it seems to be true that the time lag is 
better characterized by the centroid r ccnt of the DCF and 
the ICCF than by the peak value r pca k, namely, the time 
lag where the linear correlation coefficient has its maximum 
value r max (see e.g. iPeterson e t al. 20041. 120051 ). In both the 
DCF and the ICCF, r pca k is much less stable than r ccnt , 
but T Dea k is much less stable in the DCF than in the ICCF 
|Peterson et alj l2005h . Thus we prefer the time lag to be 
characterized by the centroid r ccnt of the ZDCF. The cen- 
troid time lag r ccnt is computed using all the points with cor- 
relation coefficients r 0.8r max in the ZDCF bumps closer 
to the zero-lag (see Fig. 3). The calculated ZDCFs between 
the radio and broad-line light curves are presented in Fig. 
3. The horizontal and vertical error bars in Fig. 3 represent 
the 68.3% confidence intervals in the time lags and the rele- 
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Table 1. Time lags between emission lines and radio emission. The sign of 



the time lag is defined 


as T ce nt — ^rad 


o — iline- Time lags are in 


units of yrs. 


Lines 


5 GHz 


8 GHz 


15 GHz 


22 GHz 


37 GHz 


Ha 
H/3 
Ho 
Ha 
H/3 
Ho 


1 9 q + 0.02 

• UZJ -0.07 
_n 1O+003 
u - 1B -0.09 
_n or+0.03 
u - 1>o -0.09 
, 79 + O.O8 
4 -' z -0.02 
4 Ql+0.07 

D - iy -0.04 


i nq+0 06 

i - My -o.oi 
-1 06+ 08 

1 1O+0.07 

J —" > — Q.Q1 

4 2n +0 ' 06 

4 - zu -0.02 
4 1Q+0.07 

4 1Q+0.06 


„ __+0.02 
z - ' ' -0.08 

-2 40+ ' 02 
z -^ u -0.07 
, r 7 +0.02 
L -°< -0.07 

n 4 n + 0.08 

q 7Q + 0.07 

,5 ' ' z -0.02 

q 7 q + 0.07 

— 0.02 


_q U +0.02 

9 O7+0.03 
ZiZ '-0.09 
i 09+O.O2 
1 - oz -0.08 
q 91+O.O7 
°- il -0.02 

O 45 + 0.07 
q 7O+0. 08 

— 0.03 


q qO+0.02 

°' ou -0.08 
q 97+O.O2 
-0.08 

-2 ni+ 002 

z - ui -0.07 
„ nfi +0.08 

q 9fi +0.09 
°'^°-0.03 
q 9 Q +0.08 
°- za -0.03 



vant correlation coefficients, respectively. The ZDCF bumps 
closer to the zero-lag have a good profile in Fig. 3. The 
measured time lags are listed in Table 1. The centroid r CC nt 
is calculated by r ccnt = E T (*) r (*)/ E r (*)> where r(i) and 
r(i) are the values of i-th data pair with r ^ 0.8r max . The 
errors of r cont are calculated by Ar* nt = {X][A T± (*) r (*) + 
r(i)Ar±(i)] E r(i) - E r(t>(i) E Ar± (<)}/E r(i)] 2 , where 
Ar ± (i) and Ar ± (i) are the relevant errors of r(i) and r(i), 
respectively. 

Our results show two possibilities that the broad-line 
variations lag or lead the radio ones (see Fig. 3). The mea- 
sured time lags are on the order of years (see Table 1). For a 
given line, the relevant time lags generally decrease as radio 
frequency increases from 5 to 37 GHz (see Table 1). The 
calculated ZDCFs between the infrared and Ha light curves 
are presented in Fig. 4 for illustration. For the ZDCFs in 
Fig. 4, there are no significant common features closer to 
the zero-lag. Also, there are no significant common features 
closer to the zero-lag for the ZDCFs of the H/3 and H7 lines 
relative to the infrared emission. On the contrary, the ZD- 
CFs in Fig. 3 have significant common features closer to 
the zero-lag. The absence of significant features tests the 
expectation in subsection 3.1. This test indicates that the 
infrared synchrotron emission does not dominate the energy 
output in the intervals between the synchrotron flares in the 
infrared-optical bands. The lag r b is the lag r ccnt measured 
here. Hereafter, T cont is equivalent to r ob and (r ob ). 



3.3 Calculations 



For 3C 273, iKaspi et al] (|2000l ) determined the Ha, H/3 
and H7 lags r r elativ e to the optical continuum, and 
iPaltani fc Tiirlerl (2005) determined their lags relative to 
the UV continuum. The optical continuum is strongly con- 
taminated by non-thermal emission, possibly related to the 
relativistic jet, and therefore it appears unsuitable for study- 
ing the lags between the ionizing co ntinuum and the lines 
i|Paltani. Courvoisier fc Walter! moj ). The Ha, H/3 and H7 
lags relative to the UV continuum are more reliable than 
those relative to the optical continuum l|Paltani fc Tiirlerl 
120051). Thus we adopt the Ha, H/3 and H7 lags determined 
bv lPaltani fc Tiirlerl (|2005h . Here, the average of rest-frame 
lags of these lines, r = 2.70 years, is adopted as a char- 
acteristic value of r. Thus the BLR has a typical size of 
Rblr = 2.70 ly. 

Simulations show that the relativistic jets can 
be driven from a region just outside the ergosphere 
of a Kerr black hole ( see e.g. iMeier. Koide fc Uchidal 
l200ll ; iKoide et all |2002| ; iKoidd |2004| ; iKomissarovl I2004T 



ISemenov. Dvadechin fc Punslvl 120041 '). In most cases, the 
bulk velocity of j et Vj is close to the escape speed 
(|Kudoh et all |l998). The escape speed is around 0.9c 
near the ergosphere of the rapidly spinning black hole 
jMeier. Koide fc Uchidal l200ll ). Most supermassive black 
holes are spinning rapidly l|Elvis. Risaliti fc Zamoraml 
2002). Thus Vj ~ 0.9c. If the disturbances in the cental en- 
gine are transported with the jet itself, Hd = fj ~ 0.9c. Thus 
we would have _R 7 ~ 9i?BLR from equation (1). 

For 3C 27 3, the jet has = 0.95c and cos 6 = 0.95 on 
100 pc scales (jDavis. Unwin. Muxlowl 119911 ). The pc-scale 
jet was constrai ned to have 8 < 1 5° and the bulk Lorentz 
factor F > 10 (|Unwin et alj 119851 ). The actual value of 9 
cannot be too small, unlike bett er aligned blazars, because 
it ha s a strong big blue bump (jShieldsl Il978l ; ICourvoisierl 
Il998h . Thus it is likely that 6 be larger than 10°. The ratio 
of the jet to counterjet flux is R = [(1 + v\ cosfl/c)/(l — 
Vi cos 9/c)] s+a for discrete moving blobs (|Lind fc B landford 



1985). R > 10 4 was observed for 3C 273 and Vj can be up to 
0.995c (see e.g. iGeorganopoulos et al.] [ 2006 ) . The observed 
spectral index a = 0.8 ( Unwin et al.lll985l ). It is obvious that 
9 < 21° and 0.9c ^ Vj ^ 0.995c are allowed by R > 10 4 . The 
disturbances are transported from the central engine down 
the jet, and then it is possible that Vd — Vj = 0.9-0. 995c 
and 9 — 12°-21°. For a given line, the relevant time lags 
generally decrease as radio frequency increases from 5 to 37 
GHz (see Table 1). We adopt the measured time lags of the 
lines relative to the 37 GHz emission. The negative lags have 



an average of t 



-2.86 years. The positive lags have an 



average of TJ b — 3.20 years. There is no zero-lag. In the 
following calculations, = 0.9-0. 995c and 9 = 12°-21° are 
adopted, and these values estimated from r D b are denoted 

by -Rradio- 

From r~ b = -2.86 years, Rblr = 2.70 ly, v d = 0.9- 
0.995c, 9 — 12°-21° and equation (7), we can obtain the 
radio emitting position i? la dio = 0.40-2.62 pc (Case B). The 
typical size of -Rblr = 2.70 ly = 0.83 pc is within the range 
of -Rradio estimated in Case B. The radio emitting regions 
in Case B are at distances of pc-scale from the central en- 
gine and are around the BLR, i.e. likely inside the BLR, co- 
located with the BLR or outside the BLR. From r+ b = 3.20 
years, Rblr = 2.70 ly, v d = 0.9-0.995c, 9 = 12°-21° and 
equation (7), we can obtain R ra dio = 9.43-62.31 pc (Case 
C). The estimated sizes are much larger than the typical 
size of i?BLR = 0.83 pc. The radio emitting regions in Case 
C are at distances of tens of pc from the central black hole 
and are far away from the BLR. 

iKovalev et ail |2009h identified the pc-scale radio core 
as a likely location for both the 7-ray and radio flares. 
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Figure 3. ZDCF between Ha and (a) 37, (ft) 22, (c) 15, (d) 8 and (e) 5 GHz; ZDCF between H/3 and (/) 37, (g) 22, (h) 15, (i) 8 and 
(J) 5 GHz; ZDCF between H7 and (k) 37, (Z) 22, (m) 15, (n) 8 and (o) 5 GHz. The x-axis is in units of yrs. 
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Figure 4. Light curves of infrared emission in (a) K, (ft) L, (c) J and (ci) H bands. ZDCF between Ho and (e) K, (/) L, (g) J and (h) 
H bands. The x-axis is in units of yrs. 



8 H. T. Liu, J. M. Bai and J. M. Wang 



Ijorstad et alj ((2001) concluded that both the radio and 7- 
ray events are originating from the same region of a relativis- 
tic jet. In 1990s, it is commonly thought that the 7-rays are 
produced in the jet, bu t closer to the central engine than the 
radio emission (see e.g. lDermer fc S chlickciscr 1993). Thus it 
is expected that Ry < i? ra dio- The constraint of i? 7 < i? ra dio 
is allowed by the recent flares of 3C 27 9 observed by Ferm i 
and in a multi-wavelength campaign |Abdoet al.ll2010dl) . 
where radio light curves from 5 to 230 GHz fail to show 
prominent variations during either the November 2008 or 
the February 2009 7-ray flares (or anytime in between). For 
Case B, Ry < 0.40-2.62 pc for v d = 0.9-0.995c and 9 = 12°- 
21°. For Case A, the zero-lag position Ry = 4.67-30.81 pc, 
which is far away from the BLR. For Case C, Ry < 9.43- 
62.31 pc. Also, Ry > 4.67-30.81 pc for the positive lags in 
Case C. Thus for Case C we have 4.67-30.81 < Ry < 9.43- 
62.31 pc for v d = 0.9-0.995c and 9 = 12° -21°. The depen- 
dence of Ry and -R ra dio on v d and 9 is presented in Fig. 5 
(see three-dimensional plots). Ry and -R ra dio increase as Vd 
increases, but decrease as 9 increases. The uncertainties of 
t>d and 9 result in the larger intervals of f? 7 and -R ra dk>- For 
better representing intervals of Ry, the sections of the three- 
dimensional plots at cos 9 — 0.95 are also plotted in Fig. 5 
(see the bottom panel). For Case B, Ry < -R ra dio = 0.44- 
1.28 pc for Vd = 0.9-0.995c and cos6» = 0.95. These Ry 
marginally satisfy Ry < -Rblr. For Case C, we have 5.15- 
15.08 < Ry < 10.40-30.45 pc. 

It is possible that there is a special point D within 
segment AG (see Fig. la). As the ionizing photons travel 
from point A to B, the disturbances travel from A to D, i.e. 
Ry = AD. Thus we have Ry/vd = Rblr/c, and then 

R, = — - — Vd- (8) 

In this case, the 7-rays will lead the lines. Combing equations 
(7) and (8), we have 

c(r ob ) 1 



-R 7 — 



(9) 



1 + 2 cos 9 

In this special case (hereafter Case D), Ry < -Rblr is ex- 
pected from equation (8). From r~ b = —2.86 years, 9 — 12°- 
21° and equation (9), we have Ry < -R ra dio = 0.77-0.81 pc. 
These estimated Ry and -Rblr = 0.83 pc satisfy Ry < -Rblr. 
This tests the correctness of -R 7 < -Rblr expected from 
equation (8). This test confirms the reliability of the time 
lags estimated by the ZDCF method. Those estimated -R ra dio 
in Case B contain these -R ra dio estimated in Case D. Thus 
Case D is a special Case B, and it is possible and reasonable. 
Combining equations (7) and (8), one can also obtain 



i^BLR. = — 



(fob) 



1 + Z Vd cos ( 



(10) 



From 



= -2.86 years, v d = 0.9-0. 995c, 9 = 12°-21° 
and equation (10), we have -Rblr. = 0.77-0.90 pc. These 
estimated values contain the typical size of Rblr = 0.83 pc. 
This confirms the reliability of the time lags estimated by 
the ZDCF method. 



4 DISCUSSION AND CONCLUSIONS 

The positions of 7-ray-emitting regions are still an open 
and controversial issue in the researches on blazars. Based 






Figure 5. Dependence of distance R on Vd and 6. From the top 
down, the first three panels correspond to Cases A, B and C, 
respectively. The x, y and z-axes are Vd in units of c, in units 
of degree and log R in units of pc, respectively. The first three 
panels represent dependence of R on v d and 9. The bottom panel 
represents dependence of R on Vd in case of cos 9 = 0.95, where 
the gridding area represents the allowed interval of iJ 7 in Case C. 
The arrows represent the upper limit of iJ 7 in Case B. 
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on the method proposed in section 2, we attempt to locate 
the emitting positions of 7-rays within the second bumps 
in the broad-band SEDs of blazars. In our previous work s 
jLiu fc Bail200d ; lLiu, Bai fc Mal200sl ; lBai. Liu fc Mal200d ). 
the internal absorption for 10 GeV-1 TeV 7-rays were used 
to constrain _R 7 , independent of how the 7-rays are pro- 
duced. Here, we try to locate R 7 , independent of the ener- 
gies of 7-rays from the SSC and EC processes. We find two 
emitting regions, the inner one at sub-pc-pc scales from the 
central black hole and the outer one around tens of pc scales. 
The outer one satisfies R 7 ^> -Rblr (Case C). The inner one 
in Case D satisfies -R 7 < -Rblr- The inner one in Case B 
mostly satisfies R 7 < -Rblr- At the same time, the inner 
one in Case B partly satisfies -R 7 > -Rblr, i.e. -Rblr < R 7 < 
2.62 pc. 

It was suggested _R 7 < -Rblr JChisellini fc Madaul 

1996). iGeorganopoulos. Kirk fc Mastichiadis (2001) argued 
R- t ;$ Rblr for powerful blazars. Tavecchio fc Mazinl (2009) 



assumed _R 7 < -Rb lr fo r the VHE 7- r ays i n 3C 279 



iLiu. Bai fc Mai l|2008l ) and iBai. Liu fc Mai (120091) suggested 



that ,R 7 is within the BLR for 3C 279. lGhisellini et all (|20ld ) 
modelled the SEDs of bright Fermi blazars, and they found 
that the position of the jet dissipation region -Rdiss is smaller 
than -Rblr for 53 out of 57 FSRQs. However, -Rdiss > R^blr 
for BL Lacs. They used -Rblr = lO 17 !/^^ cm to estimate 
-Rblr for BL Lacs and FSRQs, where Ld,45 is accretion disc 
luminosity in units of 10 45 erg s . It is appropriate to use 
-Rblr. = W 17 L^/% 5 cm to estimate -Rblr for FSRQs, but 
not for BL Lacs because this relation is derived from the 
type 1 AGNs. Thus it should be reliable that -Rdiss < R!blr 
for blazars. If the jet dissipation region is equivalent to 
the 7-ray-emitting region, -Rdiss < Rblr is equivalent to 
_Ry < -Rblr. For 3C 273, R d i ss < -Rblr l|Ghisellini et all 
|2010| ). These previous findings are consistent with our re- 
sults of < -Rblr ob tained in Cases B and D. It was 
also argued _R 7 > -Rblr |Lindfors. Valtaoia fc Turler]|2005l ; 
ISokolov fc Marscherl 1200^7 " This is marginally consistent 
with -Rblr < R 7 < 2.62 pc obtained in Case B. These 
agreements confirm the reliability of the method and as- 
sumptions. 

iBottcherl (|2008l ) suggested for 3C 279 that -R 7 > 
-Rblr for VHE 7-rays. It is recently advanced that the 
bulk of the 7-rays is produced in regions of the jet at 
distances of tens of pc from the c entral black ho l e (e.; 



Sikora. Moderski fc Madeiskil |200S| ; iMarscher et~aH 1201 .. . 
Bai fc Led l|200ll) predicted the existence of large scale syn 



chrotron X-ray jets in radio-loud AGNs, especially, the X- 
ray jets are bright on 10 kpc scales in most red blazars 
and red blazar-like radio galaxies. According to their pre- 
dictions, the large scale synchrotr on X-ray jets can p r oduce 
VHE 7-rays by the SSC process. IZhang et all (|2009l . l201fj| ) 
predicted the hot spots in lobes and the knots in jets to 
be possible GeV-TeV emitters. Fermi /L AT may resolve the 
large scale 7-ray emitters than the nuclear emitters. These 
previous findings support R-, 3> -Rblr as we obtain in 
Case C. Also, -R 7 -Rblr in Case C is not inconsi stent 
with -R 7 > -Rblr oflLindfors. Valtaoia fc Tiirler] (|2005l ) and 
ISokolov fc Marscherl (|2005l ). These confirm the reliability of 
our results. 

Our previous works are applicable to the 7-rays emitted 
from the regions in powerful blazars. where _R 7 is not much 



IBai. Liu fc Mai . 2009). The method proposed here can locate 
_R 7 in the jet. The inner emitting regions with R 7 < -Rblr. 
are likely the major contributor of the 7-rays below 10 GeV, 
for that the 7-rays above 10 GeV are subject to the photon- 
photon absorption due to the dense external soft photons 
at the inner regions (e.g. ILiu fc Bail 120061 ; ILiu. Bai fc Mai 
l200Sl ; lBai. Liu fc Mall2009l ). The outer emitting regions with 
R 1 ^> Rblr are likely the major contributor of the 7-rays 
above 10 GeV, for that these 7-rays are not subject to the 
photon-photon absorption due to the thin external soft pho- 
tons at the outer regions. Fo r these possible 7-ray emitters 
at large sc ales of kpc-Mpc (|Bai fc Led l200ll ; I Zhang et al.l 
|2009| .|2010). our works are not applicable. 

The most prominent features on VLBI images of 
jets in radio-loud A GNs are the radio co re and bright 
knots in the jet jjorstad et al.l 120071 ). iKovalev et al.l 
(2009) investigated the relation between AGN 7-ray 
emission and pc-scale radio jets. They identified the 
pc-scale radio core as a likely location for both the 7-ray 
and radio flares. A few hundreds of Schwarzschild radii, 
sub-pc-scale, is the preferr ed jet position where most 
of the dissipation occurs dGhisellini fc Tavecchid" 120091 ; 



Ghiscllini, Tavecchio fc Ghirlandal 120091 ; I Ghisellini et alJ 
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larger than -Rblr. (iLiu fc Bail 12009 ; ILiu. Bai fc Mai |200S| : 



20101 ). ISikora et all |2009l ) suggested that the blazar 
emission zone is located at pc-sc ale distances from the 
nucleus. iGhisellini fc Madaul Jl99d) s u ggest ed that R 7 is at 
sub-pc scales. Blandford fc Levinsonl ( 19951 ) also suggested 
a sub-pc 7-ray-emitting region from the central black hole. 
These previous findings support R 7 at sub-pc-pc scales. 
These i? 7 of sub-pc-pc scales are consistent with those 
-R 7 obtained in Cases B and D. These sub-pc-pc scale 
R^adio obtained in Cases B and D are also consistent with 
the previous findings of the blazar emission zone and the 
dissipation zone. These agreements confirm the reliability 
of our results. 

It is recently advanced that the bulk of the 7-rays is gen- 
erated in regions of the jet a t distances of tens of pc from 
the central black ho l e (e.g . ISikora. Moderski fc Madeiskil 
2008; IMarscher et all 12010 ). For Case C, we obtain the 
outer emitting regions of 4.67-30.81 < R 7 < 9.43-62.31 
pc and -Rradio = 9.43-62.31 pc. These outer emitting re- 
gions are comparabl e to the 7-ray-emittin g regions at dis- 
tances of tens of pc. iTavecchio et al.l l|20ld ) found evidence 
of variability on timescales of few hours from the 1.5 years 
Fermi/LAT light curves of FSRQs 3C 454.3 and PKS 1510- 
089. They concluded that significant variability on such 
short timescales disfavor the scenario in which the bulk 
of the 7-rays is produced at dista n ces of tens of pc (e.g . 
ISikora. Moderski fc Madeiskil 120081 ; IMarscher et~ai1 l201ol ). 
The previous researches show that there are two possible 7- 
ray-emitting regions, one inside or around the BLR and the 
other outside the BLR. This paper gives the same results. 
However, the method cannot discriminate between positive 
lags and negative lags on observational grounds alone (at 
least not with the current data), and the application dis- 
cussed in the paper does not distinguish between the two 
proposed scenarios. We expect this situation to change with 
future data, perhaps longer line light curves, such as 10-15 
years. The longer line light curves could give stronger con- 
straints on the coupling of the radio light curves with the 
line ones. 

For a given line, the relevant time lags generally de- 
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crease as radio frequency increases from 5 to 37 GHz. The 
trend is likely from the radiative cooling of relativistic elec- 
trons. iBai fc Led (|2003l l deduced the synchrotron time lag 
formula (see equation 9 therein). This formula can be ex- 
pressed as in the observer's frame 



1492 6 w^r("»" 2 -""' 2 )' 

where D is the "Compton dominance" (see e.g. 
iGhisellinT et alj 119981 ). B is the magnetic field strength 
in units of gauss, S is the Doppler factor, and uu and 
z^l in units of GHz are high and low fre quencies in the 
observ er's frame, respectively. For 3C 273, Ghiscllini et al.l 
(1998) obtained B = 8.9 G and 5 = 6.5. Because the 
radio light curves used to calculate the ZDCFs span 
more than 20 years and the line light curves span about 
7.5 years, it is better to derive D by using the ratio of 
synchrotron to 7-ray average luminosity. D is of the order 
of mag nitudes of 1 l|Tiirler et al l Il999n . We can obtain 
Tiag(yrs) = 12(v~ 1/2 -v~ 1/2 ) if D = 1 is adopted. The total 
cooling of both synchrotron and 7-ray emission can lead to 
Tiag(yrs) = 6(^ H 1/2 - vl 1/2 ). We calculate the ZDCFs and 
time lags between the light curves of 5, 8, 15, 22 and 37 GHz. 
The high frequency variations lead the low frequency ones. 
The measured time lags t°^ s and the relevant frequency 
differences ^ H 1,/2 — v L are presented in Fig. 6. The 
observational data are well consistent with the prediction 
of T!°g = 6(^ 1//2 — ^l 1//2 ) (see Fig. 6). This agreement 
confirms the origin of radiative cooling f or the time lags 
betwe e n the radio light curves used here. iPvatunina et al.l 
(2006, 2007) also found frequency-dependent time delays 
for strong outbursts in several other blazars. In Fig. 7, 
we compare the lags t°^ s with the differences of Aw 
between r ccnt listed in Table 1. The line of Ar ccn t = "Hag 
is consistent with the measured data points (see Fig. 7). 
This agreement confirms that the trend, i.e. the lags for a 
given line generally decrease as radio frequency increases, 
most likely results from the radiative cooling of relativistic 
electrons. 

In addition, there is another possibility that lower fre- 
quencies probe larger radii in the jet, as synchrotron self- 
absorption is important for increasingly high radius with 
decreasing radio frequency. The synchrotron self-absorption 
coefficient a v is a„ oc f~ ( - n+4 '^ 2 N B , where JV e is the elec- 
tron density and n is the electron distribution index. For a 
homogeneous blob with a radius of r, the synchrotron self- 
absorption optical depth t u is t„ = ra„ oc rv~^ n+A ^ 2 N c oc 
ru~ {n+A)/2 /r 3 = v~ {n+A)/2 /r 2 . Thus the radio frequency v 
can probe the radius r that scales as r oc ^ _ ( n + 4 ''' 4 . Hence, 
lower frequencies probe larger radii in the jet due to the syn- 
chrotron self-absorption. The higher frequencies will escape 
earlier from the blob and later the lower ones as the blob ex- 
pands. Thus the lower frequencies lag the higher ones, and 
the relevant time lags ri ag are related to frequencies. The 
difference in lags listed in Table 1 could originate from the 
synchrotron self-absorption, and it scales with frequencies 
as Ti ag oc tl — ru oc v L ( 1+n / 4 > _ ^ H (1+n / 4 \ The dependence 
of Tiag on frequencies is different from that of equation (11). 

It seems possible to infer -R ra dio based on the time lags 
of the radio synchrotron emission relative to the UV contin- 
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Figure 7. Ar ccn t vs Circles present Ho, squares H/3 and 

triangles H7. Open symbols present the negative lags and fulled 
symbols the positive lags in Table 1. Solid line is Ar ccnt = T[°g- 

uum used bv lPaltani fe Tiirleri l|2005h . This approach seems 
more direct than based on the lags of the broad lines rela- 
tive to the radio emission. We calculate the ZDCF between 
the light curves of the UV continuum and the 37 GHz emis- 
sion. There is only a little bump closer to the zero-lag for 
the ZDCF and the little bump has r max = 0.32 ± 0.06. In 
the ZDCFs between the Balmer lines and this radio emis- 
sion, the bumps used to calculate r ce nt have r max = 0.6-0.7 
that are much higher than r max = 0.32 ± 0.06. This indi- 
cates that the correlation of the UV continuum with this 
radio emission is much weaker than the Balmer lines with 
this radio emission. For 22, 15, 8 and 5 GHz, there are the 
same cases as in 37 GHz. The UV continuum is regarded as 
the ionizing continuum that drives the broad lines through 
the photoionization process. Thus it is expected that the 
correlation of the UV continuum with the radio synchrotron 
emission should be more significant than the broad lines 
with the radio emission. However, this expectation is con- 
trary to the measurements in the paper. This disagreement 
indicates that t he UV continuum is like ly not the real ioniz- 
ing continuum. iPaltani fc Tiirleri {2005) argued that the UV 
continuum is much closer t o the ionizing contin uum than the 
optical continuum used bv lKaspi et al.l j200Cf ). That is, the 
UV continuum is still not the real ionizing continuum. Thus 
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it is more reliable to derive i? ra dio from the lags of the broad 
lines relative to the radio emission than from those of the 
radio emission relative to the UV continuum. 

The relativistic shortening of variation timescales seems 
to have significant effect on the the estimates of time lags. 
The correlation between the ejection epochs of jet compo- 
nents (superluminal radio knots) and the dips in the X-ray 
emission was interpreted as accretion of the X-ray-emitting 
gas in the inner accretion disc into the central black holes 
and ejectio n of a portion of the infalling material into the 
jet (see e.g. iMarscher et al.ll2002T l. An instability in the ac- 
cretion flow causing a section of the inner disc break off. 
Part of this section is drawn into the event horizon of the 
central black hole but with considerable material and en- 
ergy ejected down the jet. The loss of this section of the 
inner disc causes a decrease in the soft X-ray flux, which is 
observed as a dip. This disturbance from the inner accre- 
tion disc to the jet is observed as the ejection of a super- 
luminal radio knot from the radio core of the jet. The dips 
in the X-ray flux represent the onset of the disturbances 
in the central engine. The knots represent the synchrotron 
emission of the transported disturbances at the sites of ejec- 
tions. Hence, the intervals between the epochs of the dips 
Aidip represent those between the disturbances in the cen- 
tral engine. The intervals between the ejection epochs of 
knots Atknot result from the disturbances in the central en- 
gine. In the radio galaxy 3C 120, it was found that there 
are 13 ejection epochs of radio knots with the c orrespond- 
ing dips in the X-ray flux (|Chateriee et al.| [2009). The ejec- 
tion epochs t kno t and the corresponding epochs of dips t d i p 
have a good correlation with tknot = tdip + 0.19. For the 
13 data pairs of tknot and tdip, At knot and Atdip between 
data pairs are equal to each other within the uncertainties, 
i.e. Atknot = Atdip- The correlation of t kno t = tdip + 0.19 
also gives Atknot = Atdip. For the ejections of knots in 3C 
120, there are the rele vant local peaks in synchrotron emis- 
sion flux from the jet l|Chateriee et al. I l2009l ;l Tavares et al.l 
|2010| ). The intervals between the epochs of peaks At pea k 
are equal to Atknot of the relevant ejections within the un- 
certainties, i.e. At pcak = Atknot- Combing At kno t = At dip 
and Atp Ca k = Atknot, we have At pea k = Atdip. The in- 
tervals Atp Ca k are measured by the light curves of beamed 
synchrotron emission from the jet. The intervals Atdip are 
measured by the light curves of un-beamed emission from 
a disc-corona system (|Chateriee et al l 120091 ). These inter- 
vals Atpeak and Atdip are generated by the same distur- 
bances from the disc to the jet. Thus it is likely that the 
relativistic effects on the time lags between them are negli- 
gible as both variations of the beamed synchrotron emission 
from the jet and the un-beamed emission from the disc- 
corona system are mainly generated by the same distur- 
bances from the disc to the jet. In fact, the DCF method 
was employed to search the time lag between the X-ray and 
the 37 GHz variations in 3C 120, and one anti-correlation 
was found wit h the X-ray leading t he radio variations by 
120 ±30 days JChateriee et alj|200a ). By the DCF method, 
ICourvoisier et al.l l| 19901 ) obtained for 3C 273 that the UV 
light curve leads the rad i o emi ssion by a few months. Based 
on this lag, ICourvoisierl (| 19981 ) showed that the radio emis- 
sion is located some 4 light years from the central source 
along the jet. Though the relativistic effects are not con- 
sidered in these DCFs used to estimate the time lags of 



the radio emission relative to the un-beamed emission of 
the UV and the X-rays, receivable results are obtained in 
these works. By analogy, it is likely that the relativistic ef- 
fects would not have significant influence on these ZDCFs 
between the radio emission and the broad-line light curves 
used in the paper. Hence, the inferred time lags from these 
ZDCFs should not be affected by the relativistic effects very 
significantly. 

For testing the correctness of the time lags obtained by 
these ZDCFs, we compare the 37 GHz light curve with the 
broad-line light curves moved horizontally and vertically (see 
Fig. 8). For the positive lags adopted, the line light curves 
are moved right by 2.8 years for Ha, 3.5 years for H/3 and 
4.0 years for H7. For the negative lags adopted, the line light 
curves are moved left by 3.3 years for Ha, 3.2 years for H/3 
and 2.0 years for H7. These moved line light curves are ba- 
sically co- varied with the radio light curve (see Fig. 8). The 
relation of At pC ak = Atdip in 3C 120 also indicates that they 
should have similar observed timescales if both variations of 
the beamed synchrotron emission from the jet and the un- 
beamed emission from the disc-corona system are mainly 
generated by the same disturbances from the disc to the jet. 
These moved times are basically consistent with those time 
lags listed in Table 1. The averages of these moved times for 
the positive and negative lags are consistent with T^ h = 3.20 



years and |t 



2.86 years obtained in these ZDCFs, re- 



spectively. The local peaks in the radio light curve basically 
have the corresponding ones in the moved line light curves. 
These indicate that the relativistic effects would not have 
significant influence on the estimates of time lags. The dis- 
turbances in the accretion flow passing through the sites 
of the central ionizing continuum could be imprinted on 
the ionizing continuum, the broad emission lines and the 
disc-jet system. In order to simulate the influence of the 
disturbances on the synchrotron emission from the jet, the 
transporting process of the disturbances from the disc to 
the jet should first be simulated by gener al relativistic mag- 
netohydrodynam ics (GRMHD, see e.g. iKomissarovl 12002 : 
lMcKinnevll2006l ). The light curves of broad emission lines 
could be simulated by the photoionization code CLOUDY 
|Ferland et al.l ll998). The GRMHD simulations might give 
the site of radio synchrotron emission, and then it will be 
easy to estimate the time lags of the radio emission relative 
to the broad lines. Also, time lags can be derived from those 
ZDCFs between the simulated light curves of radio and line 
emission. Comparing the two kinds of time lags from simula- 
tions should determine whether or not the relativistic effects 
are considered in the ZDCF analysis between one beamed 
emission and one un-beamed emission due to the same dis- 
turbances. These simulations are out of the scope of this 
paper. 

In the paper, we propose a new method to derive the 
7-ray-emitting position 7? 7 from the time lags r b of the 
7-ray emission relative to the broad lines (see Fig. 1). The 
method is also applicable to lower energy bands, such as ra- 
dio emission. depends on four parameters -Rblr, Vd, i~ob 
and 6. As r Q b = 0, r D b < and r Q b > (Cases A, B and 
C), the broad lines zero-lag, lag and lead the 7-rays, respec- 
tively. All cases are unified into equation (7). The method is 
applied to FSRQ 3C 273. Because the 7-ray light curves are 
very sparsely sampled for 3C 273, it should be unreliable to 
employ them to estimate the time lags. Fortunately, it was 
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Figure 8. The 37 GHz light curve versus the Ha, H/3 and H7 
line light curves moved horizontally and vertically, (a) the positive 
lags: moved right, (b) the negative lags: moved left. 



suggested that _R 7 < -R ra din llDermer fc Schlickeiserl 1994 



IJorstad et aLll200ll ; IKovalev et alj|2009l : 



Abdo et al 



2010dh . 



Thus R-y could be constrained by -R ra dio derived from the 
lags of the radio emission relative to the broad lines. The 
ZDCF method is used to analyze the correlations of the ra- 
dio and infrared emission with the broad lines Ha, H/3 and 
H7. The broad lines lag or lead the 5, 8, 15, 22 and 37 GHz 
emission (see Fig. 3). However, there is a lack of correla- 
tion between the infrared emission and the broad lines (see 
Fig. 4). The measured lags are on the order of years (see 
Table 1). The measured lags for a given line generally de- 
crease as radio frequency increases (see Table 1). This trend 
most likely results from the radiative cooling of relativis- 
tic electrons (see Fig. 7). The measured negative lags have 
an average of r~ nt — —2.86 years for the 37 GHz emis- 
sion relative to the broad lines. From r~ b = —2.86 years, 
i?BLR = 2.70 ly, v d = 0.9-0.995c, 9 = 12°-21° and equa- 
tion (7), we obtain i? rad io = 0.40-2.62 pc (Case B). These 
estimated -R ra dio contain the typical size of -Rblr. = 0.83 
pc, i.e. the radio emitting regions in Case B are around the 
BLR. Thus R-! < 0.40-2.62 pc for Case B and these inner 
emitting regions mostly satisfy i? 7 < -Rblr- For Case D, the 
special Case B, _R 7 < -R ra dio = 0.77-0.81 pc and i? 7 < -Rblr- 
The measured positive lags have an average of r+ nt = 3.20 
years for the 37 GHz emission relative to the broad lines. 
From r+ = 3.20 years, Eblr = 2.70 ly, v d = 0.9-0.995c, 
9 = 12°-21° and equation (7), we obtain i? radio = 9.43-62.31 
pc ^> -Rblr (Case C). Considering the zero-lag point and the 
constraint of R 1 < R ra dio, we have 4.67-30.81 < R 7 < 9.43- 
62.31 pc for Case C. These outer emitting regions satisfy 
-R 7 2> -Rblr- The dependence of R 1 and -R ra dio on Vd and 
6 and the allowed intervals of 7? 7 are presented in Fig. 5. 
-R 7 and R ra dio increase as Vd increases, but decrease as 9 
increases (see Fig. 5). The uncertainties of i>d and 8 result in 
the larger intervals of i? 7 and -R ra dio- The positions of emit- 
ting regions in Cases B and D are consistent with those a t 
sub-pc-pc scales sug ge sted bv lBlandford fc Levinsonl Jl995h , 
Ghisellini fc Madaul l|l996h . IKovalev et all (12009T I and 



Sikora et al.l l|2009h . and also with t he dissipative positions 
at sub-pc-scale proposed by others dGhisellini fc Tavecchio 



2009; 



2009; 



Ghisellini. Tavecchio & Ghirlanda 2009; Sikora et al 



Ghisellini et al.l 2010f ). The positions of emitting re- 
gions in Case C are comparable to those at distances 



of tens of pc advanced in the r ecent researches (e.g . 
ISikora. Moderski fc Madeiskil |2008l : iMarscher etafl 120101 ). 
The relative positions of the BLR and these emitting re- 
gions are also consistent with those previous findings. From 
r~ b = -2.86 years, v d = 0.9-0. 995c, 9 = 12°-21° and equa- 
tion (10), we obtain -Rblr = 0.77-0.90 pc, which contain the 
typical size of -Rblr = 0.83 pc. These agreements confirm 
the reliability of the method and assumptions. These inner 
emitting regions are likely the major contributor of the 7- 
rays below 10 GeV, and those outer ones are likely the major 
contributor of the 7-rays above 10 GeV. The method is also 
applicable to BL Lacs, in which broad lines are detectable, 
but ionizing continuum is undetectable. The BLR sizes of 
BL Lacs may be constrained by equation (10). 
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